Doping dependent quasiparticle band structure in cuprate superconductors 
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We present an exact diagonalization study of the single particle spectral function in the so-called 
t—t'—t"—J model in 2D. As a key result, we find that hole doping leads to a major reconstruction 
of the quasiparticle band structure near (ir,0): whereas for the undoped system the quasiparticle 
states near (it, 0) are deep below the top of the band near (tt/2, tt/2), hole doping shifts these states 
up to Ef, resulting in extended flat band regions close to Ef and around (n,Q). 
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Photoemission experiments on the insulating Heisen- 
berg antiferromagnet Sr2Cu02Cl2 [jjj and the under- 
doped cuprate superconductor Bi 2 Sr 2 CaCu20g + 5 
have provided unexpected results. The experimentally 
measured quasiparticle dispersion in S^CuO^C^ shows 
a band maximum at momentum (tt/2, tt/2) with a rela- 
tively isotropic effective mass, the first ionization states 
for momentum (tt, 0) are found at an energy of « O.SeV 
below the top of the band. This 'gap' at (tt, 0) then seems 
to gradually close with doping: ARPES spectra for the 
underdoped compound P^S^CaC^Os+a show that 
the parts of the band structure near (tt, 0) shift to a bind- 
ing energy between 0.2 and 0.15ey as the hole concen- 
tration increases and for the optimally doped compounds 
the quasiparticle band near (tt, 0) is within 50meV below 
the chemical potential, forming extended flat regions || . 
Taken together, and assuming that the 'effective Hamil- 
tonians' in the different compounds do not differ dras- 
tically, these data suggest a continuous but quite sub- 
stantial deformation of the quasiparticle band structure 
with doping, resulting essentially in a shift of the parts 
near (tt,0) towards lower binding energies. This scenario 
poses a formidable problem to theory in that no numeri- 
cal study of either t—J model 0] or Hubbard model 
has ever shown a behavior which is even remotely com- 
parable: for both models, the quasiparticle band is found 
immediately below Ep at half-filling and stays there more 
or less independently of doping. Here it should be noted 
that the experimentally measured binding energy of the 
quasiparticle states near (tt, 0) in the underdoped com- 
pounds H is quite large on the scale of the quasipar- 
ticle bands in the strong correlation models, which is 
~ 2J—2QQ meV. It is therefore not plausible that this 
deformation should have been simply overlooked due to 
finite size or temperature effects. 

For the undoped models the only way known so far to 
reproduce the band structure of Si^CuO^C^ within the 
framework of t— J like models is the inclusion of longer 
range hopping terms ||. To obtain a good fit to ex- 
periment requires quite substantial hopping integrals be- 
tween both second and third nearest neighbors, resulting 



in the now intensively studied t—t'—t"—J model. Antic- 
ipating that the microscopic model parameters, i.e. in 
particular t' and t", are a generic feature of the CuC>2 
plane, this seems to suggest that a description of ARPES 
experiments in terms of strong correlation models would 
require a very strong doping dependence of these param- 
eters. In this work we present numerically exact results 
for small clusters which show that on the contrary the 
t—t'—t"—J model with fixed, doping independent param- 
eters can account for the experimentally observed doping 
dependence of the quasiparticle band structure. We find 
that upon doping, the states near (tt, 0) shift towards the 
chemical potential so that there is a strongly doping de- 
pendent quasiparticle band structure and in particular 
an extended flat band region around (tt, 0) for nearly op- 
timal doping. 

We study the t—J model with hopping integrals t be- 
tween nearest neighbors, t' between 2 nd nearest (i.e. 
(l,l)-like) neighbors and t" between 3 rd nearest (i.e. 
(2, 0)-like) neighbors. We retain only the exchange in- 
tegral J between nearest neighbors. We choose t as the 
unit of energy and adopt the values J=0.3t , t'=— 0.35i 
and i"=0.25i. These parameters are within the range of 
values derived from the three band model although there 
is some variation between different estimates |1; in fact 
our main justification is the good fit to the experimental 
quasiparticle dispersion in Sr2Cu02Cl2 which they pro- 
vide (see below). We have moreover verified that our key 
result, namely the doping induced change of the quasi- 
particle dispersion, is not sensitive to this choice, and in 
particular occurs also for the case t"=0. The quantity of 
interest is the photoemission (PES) spectrum 



.4 



PES 



'k,Lu) 



1 



(H-E )- i0+ 



Cfc,„|*0>, 



where |\&o) an d denote the ground state wave func- 
tion and energy. We have evaluated this function for 
the standard 16, 18 and 20 site clusters by means of the 
Lanczos algorithm Jl^]. To begin with, Figure |l] shows 
the PES spectrum at half-filling, combining results for 16 
and 18 sites, for both the full t—t'—t"—J model and for 
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the case t'=t"= 0. We consider only the first ionization 
states, within a few J from the top of the band. As dis- 
cussed by various authors the t' and t" terms 
lift the (quasi) degeneracy of the quasiparticle states at 
(7r/2,7r/2) and (w, 0) which occurs in the t—J model, 




FIG. 1. PES spectra at half-filling for the 16 and 18 site 
cluster. Momenta marked by an asterisk refer to the 18 site 
cluster, 5-peaks are replaced by Lorentzians of width 0.03i. 
The top of the band for the 4x4 cluster is the zero of en- 
ergy, the spectra for 18 sites are shifted uniformly to produce 
a smooth dispersion. 

shifting (it, 0) to substantially more negative (on the scale 
of the total band width) binding energy. Another notable 
effect is the change in the weight of the quasiparticle 
peaks as compared to the pure t—J model: whereas the 
peak at (ir/2,ir/2) retains its original weight, the peaks 
near (7r,0), which had a comparable weight in the pure 
t—J model, become quite weak . This may well cor- 
respond to the experimental situation, where a quasipar- 
ticle peak near (n, 0) is not actually resolved E|. All in 
all, the spectra at half-filling allow for rather direct com- 
parison with the experimental results. 
We proceed to the hole doped case. Table | summa- 
rizes some properties of the two-hole ground states. For 
the 16-site cluster the absolute ground state has finite 
momentum (ir, 0). We believe that the existence of this 
state is an artefact of the special geometry of the 4x4 
cluster, which in the case t'= t"= is known jl3| to lead 
to spurious ground state degeneracies; in fact, no other 
cluster does have a finite momentum ground state for 
even hole number. The PES spectra for this state show 
a strong asymmetry between momenta along (1,0) and 
(0, 1), which also indicates that this state is not represen- 
tative for a bulk system. Finally, the lowest two hole state 
with zero momentum in the 16 site cluster (also listed in 
Table |) is nearly degenerate with the (7r, 0) state, the 
separation in energy/site being in fact much smaller than 



the variation of this quantity between different clusters. 
We therefore believe that it is reasonable to dismiss the 
(it, 0) state and consider instead the lowest zero momen- 
tum state. Next, the 20 site cluster is special in that 
already its half- filled ground state has d x z_ y 2 symmetry 
(whereas it is s for all other clusters). Adding two holes 
does not change the ground state symmetry (see Table |) , 
so that the hole pair obviously corresponds to an object 
with s-like symmetry. Thus, for all available clusters the 
symmetry of the zero momentum ground state of a hole 
pair is changed from d x 2_ y 2 for the t—J model to either 
s or p-like symmetry. More detailed analysis shows that 
the singlet s and triplet p states are nearly degenerate 
in all clusters, so that subtle geometrical features of the 
different clusters decide which one is the actual ground 
state. 

After these preliminaries, we proceed to the discussion of 
the spectral function. Figure || shows the photoemission 
spectrum for the 2 hole ground states of 16 and 18-site 
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FIG. 2. Doping dependence of PES spectra for the 16 and 
18 site cluster. The chemical potential n- for 18 sites is the 
zero of energy, the spectra for 4x4 are shifted uniformly to 
give a smooth dispersion of the 'quasiparticle peak'. 

clusters, Figure |^ presents analogouos data for 20 sites. 
The spectra for the 18 site cluster have been calculated 
for the p x -\ike ground state with S z =0. For momenta 
along (1,1), the p^-like state gives identical spectra by 
symmetry; the spectra at (2tt/3, 0) and (7r,7r/3) in prin- 
ciple could be different, but numerical evaluation shows 
that they differ only negligibly. We have chosen the zero 
of energy at the excitation energy of the first ionization 
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state, which corresponds to setting the chemical potential 
equal to the 'left side derivative' fi-=E (N) —E (N— 1). 
Assuming that the quasiparticle peaks in the cluster spec- 
tra approximately interpolate the dispersion for the ther- 
modynamic limit, this estimate for the chemical potential 
may be expected to be accurate up to the 'discretization 
error' vp\k c F — kp\ with Vp is the Fermi velocity, k c F the 
'Fermi momentum' of the cluster (i.e. the momentum of 
the first ionization state) , and kp the true Fermi momen- 
tum in the infinite system which is closest to k c F . This 
error stems from the discreteness of the fc-mesh and is 
unavoidable in any finite cluster. One may expect, how- 
ever, that it is a small fraction of the band width, and 
thus not really appreciable. Comparison with the spectra 
at half-filling then reveals quite a substantial change of 
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FIG. 3. As Figure |, but for the 20-site cluster. The zero 
of energy is /i_ . 

the band structure seen in photoemission: whereas e.g. 
the topmost PES peaks at (tv/2,0) and (tt, 0) in the 16- 
site cluster were separated by an energy of « 0.3i at half- 
filling they are now practically at the same energy. Sim- 
ilarly, the topmost PES peaks at (2tt/3, 0) and (7r, 7r/3) 
were rj OAt below the intense peak at 7r/3) in the 

18-site cluster at half-filling, but are practically at the 
same energy for the doped case. In the 20-site cluster the 
topmost PES peaks at (27r/5,7r/5) and (tt, 0) differed by 
w OAt in the undoped case but are again nearly degener- 
ate for the two hole ground state. Similarly, the PES 
peak at (7r/5,37r/5) is approximately degenerate with 
(27r/5, 7t/5) for half- filling, but is appreciably shifted up- 
wards relative to this peak in the spectra for the doped 
ground state. Irrespective of the cluster geometry and 
ground state quantum numbers the spectra for the doped 
case thus show a shift of the peaks near (jr, 0) towards 
lower binding energies so that they become (degenerate 
with) the first ionization states. This shift in energy is 
usually accompanied by an increase in peak intensity. 
This is most clearly seen in the 16 site cluster, but also 



for the 18 and 20 site clusters there is a slight increase in 
the intensity of the peaks themselves and also an upward 
shift of incoherent weight. For other momenta, such as 
(0,0) or (27r/3, 27r/3) a similar upward shift of spectral 
weight can be seen; in these cases however, the well de- 
fined peaks at half-filling disappear completely, and the 
spectra resemble incoherent continua. The doping in- 
duced deformation of the band structure may have a re- 
markable consequence, which is also illustrated in Figure 
^|: the upward shift of the quasiparticle band near (jr, 0) 
leads to the seemingly counterintuitive result that near 
Ep the decrease of electron density leads to an increase 
of electron removal weight in the fe-integrated PES spec- 
trum. It is tempting to compare this to experimental 
results on cuprate superconductors, which also show an 
increase of PES weight at Ep with increasing doping [fl4|| . 
It should be noted, that our sampling of the Brillouin 
zone is rather coarse and finite-size effects large, so that 
the fc-integrated spectra probably have little quantitative 
significance. However, the upward shift of a weakly dis- 
persive part of the band structure by 0.3 — 0Atm0.2eV 
(possibly even accompanied by an increase of its spec- 
tral weight) is a very plausible mechanism for such an 
increased PES weight near Ep. The up to now mysteri- 
ous increase of spectral weight with doping at the top of 
the band (the Fermi energy for the metals) thus is quite 
naturally explained by this model, showing that the for- 
mation of 'new states' in the charge transfer gap with 
doping is not required to explain the experiments. 
The strong doping dependence of the quasiparticle band 
structure in the t—t'—t"—J model becomes very clear in 
a comparison with the 'pure' t—J model, shown in Figure 
[|. When compared with Figure ^, this figure shows the 
remarkable 'rigidity' of the band structure: the quasipar- 
ticle peaks in the interior of the antiferromagnetic Bril- 
louin zone (AFBZ) persist with practically unchanged 
spectral weight and dispersion, peaks outside the AFBZ 
loose most of their weight, but persist as 'shadow band'. 
For momenta near (tt, 0) the PES peaks cross to the 
inverse photoemission spectrum with almost unchanged 
weight (the fact that this happens near (it, 0) rather than 
(tt/2, 7t/2) is due to the d x 2_ y 2 symmetry of the two-hole 
ground state M). In cases where there is an apprecia- 
ble change of the peak positions, the PES peaks usually 
shift away from the chemical potential as compared to 
the strong peaks right at Ep; moreover, the weight of the 
peaks usually decreases. The general trends with doping 
therefore are opposite to that of the t—t'—t"—J model. 
Moreover, the spectra for the pure t—J model on the av- 
erage look 'more coherent', i.e. the decay of sharp peaks 
into incoherent continua is much less pronounced as for 
the t—t'—t"—J model. As far as a 'dispersion' can be 
assigned, the band for the t—t'—t"—J model essentially 
approaches that of the ordinary t—J model, i.e. t' and t" 
seem to be renormalized towards zero with doping. This 
suggests a loss of 'spin coherence' as the source of the 
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doping induced band deformation: in a Neel state, 
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FIG. 4. Doping dependence of PES spectra for the 16 and 
18 site cluster for the t—J model. 

transferring a hole to a second- or third-nearest neighbor 
gives precisely the same state as if the hole had been cre- 
ated there in the first place. The t' and t" terms therefore 
contribute directly to the coherent motion of the hole, 
and consequently have a strong impact on the dispersion. 
On the other hand, when the spin correlation length is 
comparable or shorter than the 'range' of the t' and t" 
terms, a spin transferred by these may no longer 'fit in' 
to its new environment, so that the hopping process with 
a high probability produces spin excitations and conse- 
quently becomes incoherent. In other words, rather than 
contributing to the coherent hole motion, the t 1 and t" 
terms now enhance the emission of spin excitations. This 
would also explain the generally 'more incoherent' char- 
acter of the PES spectra in the t—t'—t"—J model. By 
similar arguments the introduction of frustration at half- 
filling should also lead to a deformation of the quasipar- 
ticle band with long distance hopping, as indeed seems 
to be the case Jl5[ . 

In summary, we have studied the doping dependence of 
the photoemission spectrum for the t—t'—t"—J model. 
Using parameter values which give a realistic fit to the ex- 
perimental quasiparticle band structure in the undoped 
Heisenberg antiferromagnet Sr-iC^CliOi we found a 
substantial doping dependence of the topmost photoe- 
mission peaks: upon doping, the quasiparticle band near 
(•7T, 0) moves to lower binding energy, and becomes degen- 
erate with the first ionization states. This resuls in par- 
ticular in an increase of fe-integrated spectral weight near 
Ep with doping. While exact diagonalization does not 
allow for a sufficiently smooth variation of the hole con- 
centration to map out the detailed doping dependence, 



the rather substantial deformation of the quasiparticle 
band without changing any microscopical parameter of 
the model suggests a unifying point of view for both 
Sr2Cu02Ci2 and cuprate superconductors. As the mech- 
anism for the doping induced deformation we propose 
the interplay between longer range hopping and decreas- 
ing spin correlation length. As Neel order is lost over 
the range of the t' and t" terms, the effect of the latter 
changes from coherent hole propagation to emission of 
spin excitations, leading to the deformation of the band 
structure. Taking this further, one might expect that 
a degradation of the spin correlations due to increasing 
temperature has an analogous effect, i.e. an overall 'flat- 
tening' of the band structure. The resulting increase of 
the band mass with temperature then could provide a 
very simple explanation for much of the 'spin gap' be- 
haviour observed e.g. in the spin susceptibility and en- 
tropy of cuprate superconductors [jltjf . 
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TABLE I. Two hole ground state properties for various 
cluster sizes: momentum feo, point group symmetry R, total 
spin Stot and energy/site Eo/N. 
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